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We generalize the Swiss-cheese cosmologies so as to include nonzero linear momenta of the asso- 
ciated boundary surfaces. The evolution of mass scales in these generalized cosmologies is studied 
for a variety of models for the background without having to specify any details within the local 
inhomogeneities. We find that the final effective gravitational mass and size of the evolving inho- 
mogeneities depends on their linear momenta but these properties are essentially unaffected by the 
details of the background model. 
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I. INTRODUCTION 



The Swiss-cheese models give us (noninteracting) in- 
homogeneities in a cosmological setting that are exact 
solutions to Einstein's equations. As a result, the mod- 
els have become a standard construction [l| and are very 
widely studied [2[. The classical Einstein-Straus vacuole 
model, which requires a comoving boundary surface, is 
unstable. Many subsequent Swiss-cheese models have 
also assumed that the associated boundary surfaces re- 
main comoving, but it is well known that this need not 
be the case. Some general studies of this issue go back 
many years Q • Some specific examples of noncomoving 
boundary surfaces include the study of Vaidya-type in- 
homogeneities [J] and the evolution of "density waves" 
in the Lemaitre-Tolman model Q. Here we examine 
the role of the linear momentum of a boundary surface 
in a Robertson- Walker background. Within the context 
of these generalized models, we study the evolution of 
the effective gravitational mass and size of the inhomo- 
geneities for a variety of well-known background models. 
We find that whereas the momentum plays a central role, 
the details of the background model are relatively unim- 
portant. 



Junction Conditions 



II. MODEL CONSTRUCTION 



The theory of hypersurfaces in spacetime is well estab- 
lished (see, for example, [y]) and we do not reproduce all 
the necessary machinery here. Rather, we go directly to 
the essential ingredients of the model. The model consists 
of randomly distributed nonintersecting spherical bound- 
ary surfaces E in a Robertson- Walker background with 
each particle on S executing radial timelike geodesic mo- 
tion. 



To establish notation, let us write the metrics of the 
spacetimes V± in the form 



ds 



± 



dsl ± (x 1 ± ,x 2 ± ) + R 2 (x] t 



c 2 ± )dfl 2 



(1) 



where the signature of the two-surfaces r± is zero and 
dfl 2 is the metric of a unit two-sphere, which we write in 
the usual form d9 2 + sin 2 (9)d(f> 2 . We need consider only 
one boundary surface. The metric on £ can be written 
in the form 



ds| = R 2 (r)dti 2 



dT 2 



(2) 



where r is the proper time on S. Since £ is, by assump- 
tion, geodesic, there is only one nonvanishing indepen- 
dent component of the extrinsic curvature K a p. This 
can be written in the form (e.g. 0) 



K 2 



R 2 R 2 



1 - 



2M-_ 



(3) 



where ■ = d/dr and Ai± are the effective gravitational 
masses of V± . The invariant properties of M. were first 
explored by Hernandez and Misner Q who wrote the 
function in the form 



R 3 
2 6 



(4) 



where 1Z is the Riemann tensor of V. Sec also 
for further discussion [l3j . From ([3]), and the continu- 
ity conditions, we arrive at the central condition of the 
model: 



M- = M 



(5) 



a statement which is independent of the coordinates 
(x±,x ± ). 



B. Background Geodesies 
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The Robertson- Walker geometry (excluding the Ein- 
stein static subcase) has a Killing algebra of dimension 
6 (3 translations and 3 rotations). Since £ is in radial 
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motion, we are interested in the constants of motion gen- 
erated by translational invariance. Write the background 
in the form 



and 



deli = a 2 {t) 



dr 2 



+ r W )-dt 2 , (6) 



where e(r) = vT— kr 2 with k = ±1, 0. The translational 
Killing vectors are then given by 



iV = < [ cos(8)d r - ^-d e ) . 



e 2 Q = e 



and 



(7) 



sin(0) sin(0)a r + - (cos(e) sin(0)5 e + ^§-d, 
r \ sm(t') 

(8) 



sin(0) 
' sin(0) °* 
(9) 

= fd r + tdt 



(10) 

T we have 



sin(0) cos((/))d r H — cos(0) cos(4>)de 



Consider a radial geodesic with tangent u a 
and define the constants C n = £"u Q and T> 2 ; 
follows directly from Q-© that 

.o 2? 2 e 2 



so that from the timelike condition u a u a = 

•2 ^ 

i 2 = 1 + —f. 



(11) 



Equations (fTUT) and (TTT1) should be well known [lj| ■ Our 
purpose here is to explain the physical meaning of T>. It 
is the total linear momentum of S [15| . 



G a p n a n p + A = 0. 



(15) 



From ([6]) and (fT4|) we find r = and from ([6]) and (fT5|) 
we find 87rp = where p is the comoving isotropic pres- 
sure. As a result, in the Einstein-Straus model, S must 
be exactly comoving and the background must be exactly 
dust. If either of these conditions do not hold then ((5]) is 
necessarily violated. Without resorting to surface layers, 
we can take the view that the culprit is the assumption 
that V_ is exactly vacuum. That is, (fT3|) is too strong 
a condition to impose on V_ . In other Swiss-cheese cos- 
mologies, (|13p is not imposed on V_. A widely used ex- 
ample is to assume that V_ is dust. However, typically 
E is taken as comoving a priori and so we arrive back at 
(Q~4|) and (fTS"j) which, we would like to emphasize, are not 
necessary conditions for a boundary surface. 



D. Generalized Model 

In this paper the only condition imposed on V_ is ([5]) 
and we concentrate on the evolution of E in V+ which, for 
the sake of clarity, we assume is spatially flat. It follows 
from (HO]) and (HT) that [TJ 




V 2 (l 



dx 



V 2 (l 



H(x) 



(16) 



Swiss-cheese 



where we have set the initial conditions by r E (t = 0) = 0; 
the Hubble function H is given, as usual, by -4f; and, 
without loss in generality, we have set a (today) — 1. We 
have chosen as an independent variable z — 1/a — 1 and 
so for universes with a big bang that do not recollapse 
— 1 < z < oo. 



The standard Swiss-cheese inhomogeneous cosmology 
takes T> = 0. Of these models, the Einstein-Straus case, 
which sets V_ as vacuum, is the most well known. The 
model has seen very wide application [l| , but it is known 
to be unstable: to aspherical perturbations [lH and to 
perturbations in the condition ([5]) which could lead to 
the development of surface layers [l?} ■ A clear way to 
see the instability in this model at a primitive level is to 
look at the necessary conditions for a boundary surface: 

[G>X] = [G%n a nP] = 0, (12) 

where is the Einstein tensor, = (* + - *_)|s, 
and u a and n a are the tangent and normal vectors to E 
respectively. By assumption, V_ satisfies 

G% + A5p = , (13) 

where A is the cosmological constant. As a result, the 
following two conditions must be satisfied in V+ if E con- 
stitutes a boundary surface: 

GZriavP = 0, (14) 



III. BACKGROUND MODELS 

In this work we explore mass scales associated with var- 
ious models for the background specified by H{z). The 
inhomogeneities considered here are governed by one free 
parameter, the total linear momentum T>. The first task 
at hand then is to establish a reasonable range of values 
for this parameter. 

A. The Range in V 

To establish a reasonable range in T> we evaluate 
r E (0) from (fT6|) for three standard models: ACDM, 
ACDM+noninteracting radiation (see the Appendix) and 
noninteracting matter and radiation without A. The re- 
sults are shown in FigJT] The top curve is included so as 
to show the effect of ignoring the background radiation, 
and the two bottom curves are included so as to show 
the effect of ignoring A and the effect of changing the 
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integration upper limit to a redshift where radiation can 
be neglected, here z — 1500. The middle curve shows 
the size of the hole for the ACDM+noninteracting ra- 
diation at earlier time r s (100). The wide difference in 
size r s between the top and bottom curves shows the ef- 
fect of adding radiation on the evolution of the surfaces 
for a given value of T>. From the results of Fig. [TJ wc 
have chosen to include radiation (see the Appendix) to all 
background models H{z). Throughout the paper, we de- 
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FIG. 1: The abscissa is the dimensionless fraction T>/c x 
10 _J and the ordinate gives r E (0) in units of Mpc. The 
top curve is the ACDM model. The middle curve is 
ACDM+noninteracting radiation (with the curve just below 
giving r s (100) for comparison). The second to bottom curve 
is the noninteracting matter and radiation without A, and 
the bottom curve is the ACDM with integration upper limit 
going to z = 1500. We have adopted the values il mo = 0.27, 
h = Ho/100 = 0.72 and T = 2.73. 



fine quantities associated with a flat "standard" ACDM 
model as " ~ " when using the following values: for the 
matter energy density, Q mo — 0.27, radiation energy den- 
sity (including relativistic neutrinos), £Ir ~ 8x 10~ 5 and 
cosmological constant SI a, regardless of the value of our 
free parameter T>. Using the ACDM background, Fig. [5] 
shows the evolution of r E with total linear momentum 
T> in the order of ~ 10~ 5 . These inhomogeneities begin 
with f s (oo) = and grow mostly during the radiation 
epoch. We are interested in ranges of the parameter D 
where surfaces r s reach up to ~ 100 Mpc. From Figs. Q] 
and [21 this correspond to a range 10 -5 <P<5x 10~ 5 . 
In subsequent figures we use T> — 5xl0~ 5 unless specified 
otherwise. 



B. Flat standard cosmological model (ACDM) 

In the flat ACDM model, dark energy is the cosmolog- 
ical constant, and as usual we have 



H(z) 



n m (i + zf + n A + n K (i + z) 4 



(17) 



where the Jl refer to current values. Note that we have ex- 
plicitly included radiation density £1^ (see the Appendix) 
and it is included implicitly in all the following back- 
ground models. The influence of VL m on r E is shown in 
Fig. H The WMAP5 [3 value of fl mo = 0.279 ± 0.015 
has limited effect on the final size compared to the value 
of the linear momentum in Fig. [2] To compare the dif- 
ferent backgrounds, we can plot the mass ratio between 
the background model M.^ and the ACDM background 
M s . From (gj) it follows that 



(18) 



Throughout, for (fla)) . we use the same value for the 
linear momentum {T> — T>). For the ACDM model, 
the mass ratio (|18l) is shown in Fig. @] where A4 is 
determined by fl mo — 0.33 and different values of 
T> are shown. In Fig. [5] wc plot the mass ratio (|T5|) 
for the ACDM model for the same values of f2 mo as Fig. [3] 



ACDM 
D = 9xl0" 5 
D = 7xl0" 5 
D = 5xl0" 5 
D = 3xl0" 6 
D=10- 5 




FIG. 2: Variation in the evolution of r E with redshift horn 
variations in T> for a ACDM background. We have used 
h = 0.72 and To = 2.73, and like all "~" quantities we have 



used fl„ 



0.27 and Q Ro 



x 10" 



These surfaces start 



with r s (z — 00) — 0. The value of the linear momentum T> 
determines the final size of these inhomogeneities. 
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Q =0.27 
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FIG. 3: Evolution of r s with redshift for some values of the 
matter density fl m . We have used a ACDM background with 
a total linear momentum of V = 5 x 10~ 5 and we have extrap- 
olated the sizes of these inhomogeneities to the infinite future 
(z = —1). The dotted line shows the size at the present red- 
shift (z = 0). We have used values of fi ra consistent with 
the 99.9% confidence interval of [2(| to show the effect of 
Q,m - For comparison, the latest WMAP5 results give 
n mn = 0.279 ± 0.015. 



D. Flat dark energy with variable equation of state 

(VES) 

Allowing the equation of state to vary with time 
(VES), and using the parameterization [2l| 



w{z) — w + w a 



1 + z' 



we have 



(^) 2 =a " (1 + + nde{1 + z) Hl+W0+Wa) 



(21) 



(22) 



x e 



-3u/ a (*/(!+«)) 



The situation is examined in Fig|5] Whereas the mass 
ratio decreases with increasing w a , there is no evidence 
that w Q + w a > [22 1 . 



E. Flat DGP models 

The flat Dvali-Gabadadzc-Porrati DGP model [H| (see 
also (24[) is a one parameter model from brane theory 
where ft r — l/(4r 2 i?Q) is the dimensional parameter de- 
termined by the scale length r which governs the transi- 
tion from 4D to 5D behavior. For this model, the Hubble 
parameter is given by 



(23) 



Define 

AM = M B /M S -1. (19) 

For r E ~ f E then AAi ~ 3Ar, where Ar = (r E /f E ) — 1 
is the size difference for a fixed V. This can be gleaned 
from Figs. ME 



C. Flat dark energy with constant equation of state 

This model (which we designate by CES) uses a con- 
stant arbitrary value for the equation of state parameter 
w = p/p so that the Hubble function becomes 




= (l m (l + zf + n de (l + zf( 1+w \ (20) 



In Fig(5]we show the ratio Ai^/Ai^ for w = —0.7 and 
w = —1.3 corresponding to a 99.9% confidence level [2fj| . 
The mass difference AAi between the two models is less 
than 2% at maximum. The difference would be negligible 
using WMAP5 19] where w is constrained to — 0.97±0.06 
as can be seen in Fig[7] where we have plotted the mass 
ratio for different values of w with Q m = 0.27. 
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FIG. 4: Evolution of ratio M^/M^ with redshift for some 
values of T> = T> where M is given by the ACDM with Q m — 
0.33 and M is our standard ACDM model. 
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w=-0.8 
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FIG. 5: Evolution of ratio M s /M s for ACDM with V = FIG. 7: Ratio M s /M s for V = 5 x 10" 5 for different values 
5 x 10 -J for some values of f2 m . of w with Q m ~ 0.27 for the constant equation of state model. 
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Z 

FIG. 6: The same as Fig. [4] where A4 now has the constant 
equation of state model (CES) for the background. 

FIG. 8: The same as Fig. [7] with the variable equation of state 
(VES) model as background for M.. Only the region w a < 1 

is of physical interest. 

where f2 m = 1 — The 99.9% confidence level of 

20] was used for the values of O r in order to show the 
mass ratio (fT5|) for the DGP model as given in Fig. |H1 
Subsequent models showing the same type of behavior 

as Fig. IH1 over the redshift range — 1 < z < 20 will be summarize in Table fl] 
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TABLE I: Model properties 



Model 






1 st parameter 


2 nd parameter 


M s /M s (*) 
a = -1 


^ = 20 


AFF 


0.27 


i - n 


a = 0.02333 




0.696 


0.692 




0.27 


1 — 


n — —0 01 667 




1.245 


1.248 


CHA 






4 = 0.7 


7 = 0.25 


0.778 


0.783 








A = 0.75 


7 = -0.1 


1.171 


1.166 


IDE 




0.73 


£ = 2.3 


to = — 1 


1.306 


1.297 






0.73 


£ = 3.3 


TO = —1 


0.9192 


0.9211 


NADE 






n = 2.5 




0.859 [25J 


0.864 








n = 3.1 




1.1125 [25] 


1.1118 




FIG. 9: Same as Fig. [7] M is given for the Dvali-Gabadadze- 
Porrati (DGP) model. 



FIG. 10: Ratio M s /M s evaluated at z = —1 as a function of 
£ for w = — 1, T> = 5 X 10 _J and O^e = 0.73 for the interacting 
dark energy (IDE) model. 



F. Flat interacting dark energy 

It is natural to consider the coupling between dark 
energy and matter and there are many explicit coupling 
procedures considered in the literature. Here we use the 
parameterization of 26] to write 

fH(z)\ 2 , N ,[ / 1 \]- 3 (-/«) 

(24) 

Figure [TU] shows the mass ratio (|T5)) for the interacting 
dark energy (IDE) model with w = — 1, Sl^ = 0.73 and 
2? = 5xl0~ 5 atz = — 1 over a range of £ consistent with 
[26( | . The mass ratio (fT5)l does not change significantly 
over the redshift range — 1 < z < 20 as can be seen in 
Table D 



G. Cardassian models 

The modified polytropic Cardassian models (hereafter 
CAR) are three parameter models that modify the Fried- 
mann equation in a flat, matter-dominated universe in 
order to allow acceleration. The Hubble function is given 
by (see [HUlI) 

(&j 2 = n m (i+zf (1 + W - i)(i + z)^ n ^) 1/c 

(25) 

The evolution of the mass ratio (fl~8f with redshift for the 
Cardassian models is shown in Fig. [TT] for VL m — 0.27 
and T> = 5 x 10~ 5 and values of the parameters n and q 
surrounding the 99.9% confidence level of [20I [27|. 
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FIG. 11: Same as Fig.[7]but now M is given for the cardassian 
(CAR) model. 

H. Flat Chaplygin gas 

A dark fluid that combines dark matter and dark en- 
ergy, where the equation of state is p — —A/p 1 [28], gives 
the generalized Chaplygin gas equation (CHA) 

(igil) 2 .^ + (1 - A)a + ^»)""->, (26) 

Since the Chaplygin gas is a two parameter model we 
choose two extreme values for A for Fig. [T^] in which we 
plot the mass-ratio {JXSJ) at z = —1 as a function of 7. 
The values chosen here follow the 99.9% confidence level 
of 20]. As with the IDE model, the mass ratio for the 
Chaplygin gas model is nearly constant over the range 
-1< z < 20. See Table Q] 

I. Flat affine equation of state 

The assumption that dark energy and dark matter are 
a single dark component that can be modeled by the 
affine equation of st ate p = po + ap (AFF) gives rise to 
the Hubble function [29| 

(^) 2 = n m (i + zf( 1+a) + n A , (27) 

where Q m = (po — pa)/ Pc- The ratio (fl"8| for the AFF 
model is shown in Fig. [T3]at z = — 1 as a function of a 
using the 99% confidence level of [29( for a. Complemen- 
tary information is available in Table HI 



FIG. 12: Same as Fig. HQ] but with A = 0.67 and 0.78 as 
a function of the parameter 7 for the Chaplygin (CHA) gas 
model. 




FIG. 13: Same as Fis.HOlbut for the affine (AFF) parameter 
model as a function of a. 



J. New agegraphic dark energy (NADE) 

A single-parameter model where the energy density of 
quantum fluctuations of Minkowski spacetime p q is in- 
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eluded in the Hubble equation gives [3C 

2 



where from (IT6 



H(z) 



m mo (i + z) 3 
i-n q {z) 



(28) 



where n mo = 1 — f2 g (0) and the evolution of fl q (z) is 
given by 



dz 



-n q {i-Q g ) 



3(1 + Z)- 1 



'fin 



(29) 



The mass ratio ([T5)l for this model is presented in Fig. [HI 
with values of n consistent with the likelihood values in 
301. 




FIG. 14: Same as Fig. [7] M is given for the new agegraphic 
dark energy (NADE) model 



IV. DISCUSSION 



' V 2 {l + z) 2 1 
l + T> 2 (l + z) 2 Wz) 



(32) 



and equivalently for the ACDM model (with a bar). 
Therefore the mass ratio is extremal for 



H 
H 



My 



(33) 



An extremum at z = would imply that r E = r s and 
that M s = M. s . Equation (j3"Tj) can also be expressed as 



= 3<5 3 



V 2 



1 



l + V 2 r s (0)H 



{So - 1} , (34) 



where Sq = (Al s (0)/Al s (0)) and r E (0) is obtained from 
(fT6| for the model investigated. 

In Eq. (fTB| . the influence of 2? in the first factor of the 
integral converges to 1 for T> 2 (l + z) 2 3> 1 or 1 + z > 1/2? 
and converges to 2?(1 + z) when 1 + z <C 1/2?. In effect 
then since 2? ~ 1CP 5 , T> is unimportant at early times for 
7V and 7V, . 



V. CONCLUSION 

We have introduced a new generalized Swiss-cheese 
model which does not assume a priori that the asso- 
ciated boundary surfaces are comoving. In order to 
quantify evolving inhomogeneities, we have considered 
geodesic boundaries characterized by their linear momen- 
tum T>. For the size of the inhomogeneities we are inter- 
ested in, the physical values of V/c are ~ 10~ 5 . For a 
given linear momentum, we have found that the inhomo- 
geneities grow almost independently of the background 
model (with the inclusion of the radiation density param- 
eter fifl ). As shown in FigFi] these inhomogeneities are 
almost at their full size by the decoupling (z ~ 1100). For 
a redshift of z < 2, corresponding to high redshift super- 
novae, the inhomogeneities considered here are growing 
very slowly as is shown in see Figj4] 



A property of the mass ratio (fT8| seen in all relevant 
figures is simply 

H(z)^H(z)^M s ^M s . (30) 

As can be seen in Fig. [6j [7] and [Til t ne mass ratio (fT8|) 
appears to have an extremum at z = 0. (This feature 
is present in all plots of the mass ratio as a function of 
z but is not visible in all figures because of the chosen 
scale.) In fact, the extremum can not happen at z = 0. 
From (fTrJ)) and (fTS]) , and defining ' = d/dz, we find 

(^)' = 3 | Rr " S " f ^ } ' (31) 



ACKNOWLEDGMENTS 

CG is supported by the Ontario Graduate Scholarship 
in Science and Technology and the Fonds quebecois de 
recherche sur la nature et les technologies. KL is sup- 
ported by a grant from the Natural Sciences and Engi- 
neering Research Council of Canada. Portions of this 
work were made possible by use of GRTensorIl\3^ . 

APPENDIX: RADIATION 

The integration in Eq. (fTB)) requires a high redshift con- 
tribution from radiation. See, for example, p^l[3l|. This 
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was added to all the Hubble background models H(z) 
as noninteracting species oc a~ 4 . The energy density for 
radiation pn is 



PR = a B T ( 



CMB 



1 + U^L.)\ 
8 V Tcmb / 



(35) 



where as = 4er/c in which a is the Stefan-Boltzmann's 
constant, Tcmb = 2.725-ftf is the temperature of the cos- 
mic microwave background (CMB), T v = (4/11) 1 / 3 Tcmb 
is the temperature parameter for the relativistic neutri- 
nos after the annihilation of electron-positron pairs and 
JV„ = 3 is the standard number of neutrino families. The 



first right hand side term is the contribution from the 
CMB photons and the second term is the contribution 
from relativistic neutrinos. The energy density of radia- 
tion is therefore: 



PRo 
Pea 



X 10" 



(36) 



where p c is the critical energy density; the subscript is 
used to specify present values (z = 0). Note that before 
the electron-positron pairs annihilated (z ~ 10 10 ) the 
temperature of the neutrinos and CMB radiation was in 
equilibrium T v — Tcmb- This contribution in Eq. ([33)) 
was negligible in our integration of (|16[) and was ignored. 
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